tH 


CALCUUTIONS  OF  SOME  UNSTEADY 
^ TRANSONIC  FLOWS  ABOUT  THE  NACA 
^ 64A006  AND  64A010  AIRFOILS 


GENERAL  DYNAMICS  CONVAIR  DIVISION 
KEARNY  MESA  PLANT 
5001  KEARNY  VILLA  ROAD 
SAN  DIEGO.  CALIFORNIA  92138 


JULY  1977 


TECHNICAL  REPORT  AFFDL-TR-77-46 

FINAL  REPORT  FOR  PERIOD  NOVEMBER  1975  - APRIL  1977 


Approved  for  public  rcletMc;  distribution  unlimited 


Prepared  for 

&Sair  force  flight  dynamics  laboratory 


Air  Force  Wright  Aeronautical  Laboratories 
Air  Force  Systems  Command 
Wright'Patterson  Air  Force  Base,  Ohio  45433 


D D C 

“lEiBEnnjiE 


NOV  3 1977 


1SEDU15 

B. 


r 


Unclassified 


SeCURITV  CLASSIFICATION  OF  THIS  PAGE  (Whin  Dmtm  Enffd) 


'/o 


^7)  REPORT  DOCUMEhITATION  PAGE 

».  WEPQRy^BER  / |2.  GOV 


AFFDLiyTR-77-46  j 


r 


^TgOVT  ACCESSION  NO. 


4.  TITLE _ 

Calculations  of  Some  Unsteady  Transonic  Flows  About/ 
the  NACA  64A(i^a6  and  64AP10  Airfoils  , p''  


7.  authorc*; 


DrAR.  J. 


a gnus  J 


(/^ 


' 


*.  PERFORMING  ORGANIZjyriON  NAME  ANO  ADDRESS 

General  DynamicyConvair  PsAFtoiow- 
P.  O.  Box  80847^ 

San  Diego,  California  92138 


{£{ 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Air  Force  Flight  Dynamics  Laboratory,  FBR 
Air  Force  Wright  Aeronautical  Laboratories 
Air  Force  Systems  Command 
Wright- Patterson  Air  Force  Base.  Ohio  45433 


14.  MONITORING  AGENCY  NAME  A AODRESSCfl  dll/Aranl  /root  Conirollint  Oilier) 


A 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


CIPiENT*S  catalog  NUMBER 


.T  6 PERIOD  COVERED 
rgt  for  ppriori .... 

77^ 


75— -Aprl 


S-  PERFORMING  ORG.  REPORT  NUM 


0.  CONTRACT  OR  GRANT  NUMBER<i> 


F33615-76 


-C-3^18 


10.  program  ELEMENT,  PROJECT.  TASK 
AREA  0 WORK  UNIT  NUMBERS 


61102F  A / , 


} July,-I«77 


TJ.  NUMBER  OF  PAGES 

88  + ix 


IS.  security  class.  (oI  Ihir  raportj 

Unclassified 


ISa.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ol  llllm  Raparu 


Approved  for  public  release;  distribution  unlimited. 


D D C 

n]Ef?0F?nn  .nrp| 

Uj  NOV  3 1977 


17.  distribution  statement  (ol  Iho  aSalracI  aniarad  In  Block  30,  II  dllloronl  (root  Raporr; 


nnb'Eij'Li  A 

B 


IB.  VUPPLEMENTARV  NOTES 


10.  KEY  WORDS  fConriiMM  on  rovoroo  oJdo  1/ nocooaory  and  Identity  by  block  number) 

Transonic  Airfoil  Flow 
Unsteady  Flow 
Finite  Difference 


Unsteady  transonic  flows  over  airfoils  were  calculated  using  a program  based  on 
the  unsteady  Euler  equations.  The  approximate  numerical  solutions  were  obtained 
using  an  explicit  (Lax-Wendroff)  difference  scheme. 

Calculations  were  made  for  the  64A006  airfoil  at  zero  angle-of-attack  in  flows  with 
Mach  number  0.  822,  0. 854  and  0. 875.  Three  unsteady  flows  caused  by  a quarter-chord 
flap  oscillating  with  an  amplitude  on  the  order  of  one  degree  at  specific  reduced  ^ 


DD  , 


^ORM 
JAN  73 


EDITION  OF  I NOV  SS  IS  OBSOLETE 


Unclassified 


security  CLASSIFICATION  OF  THIS  PAGE  flWnn  Duo  Enitrod) 


/'/  7 


/ ~o 


iiMBilhiiH 


7 

_ - — p frequencies  (k  = wC/U„)  near  0.4  were  analyzed.  The  results  were  compared  with 


other  available  calculations  and  the  experimental  data  of  Tljdeman. 

Responses  for  an  airfoil  in  a ventilated-wall  wind  tunnel  rather  than  in  a free- 
stream  were  also  calculated.  Exploratory  stud^  of  the  changes  to  be  expected  with 
a weakened  shock  (an  attempt  to  simulate  a main  effect  of  the  interaction  of  the  shock 
with  the  airfoil  boundary  layer)  was  done  for  one  case  at  Mach  number  0. 875. 

Oscillatory  flows  over  the  64A010  airfoil  at  zero  angle-of-attack  in  a Mach  0. 80 
stream  were  also  calculated.  Plunging  motion  at  reduced  frequency  0. 4 and  pitching 
motions  at  reduced  frequencies  0.4  and  0.5  were  calculated.-^ 


FOREWORD 


This  report  was  prepared  by  the  Convair  Division  of  General  Dynamics, 
San  Diego,  California,  for  the  Analysis  and  Optimization  Branch  of  the  Air 
Force  Flight  Dynamics  Laboratory  (AFFDL),  Wright- Patterson  Air  Force 
Base,  Ohio.  The  research  was  conducted  under  Contract  F33615-76-C-3018, 
"Time  Dependent  Finite-Difference  Solutions  for  Transonic  Flow.  " 

Dr.  James  J.  Olsen,  FBR,  of  AFFDL  was  the  Air  Force  task  engineer. 

At  Convair,  Dr.  H.  Yoshihara  was  the  principal  investigator  during  the 
period  from  November  1975  through  October  1976.  The  preparation  of  com- 
puter programs  and  the  calculations  were  carried  out  by  Dr.  R.  J.  Magnus. 

This  report  was  submitted  in  April  1977  to  cover  research  performed 
from  November  1975  through  April  1977. 

The  support  of  this  research  by  Harvard  Lomax,  Charles  Coe,  and 
Sanford  Davis  of  the  NASA  Ames  Research  Center  and  the  furnishing  of 
calculations  and  data  for  comparisons  by  William  Ballhaus,  Peter  Goorjian, 
Joseph  Steger,  and  Louis  Stivers  of  NASA  Ames  are  gratefully  acknowledged. 


ill 


II 


III 


IV 


yhece^lxy^^  “r^A  ^ar/r/f  - 


TABLE  OF  CONTENTS 


SECTION 

I INTRODUCTION 


METHOD  OF  COMPUTATION 


PAGE 

1 


1. 

Equations 

3 

2. 

Grid  Arrangement 

4 

3. 

Difference  Scheme 

5 

4. 

S5mchroniz  ation 

6 

5. 

Boundary  Conditions 

6 

6. 

Computational  Notes 

8 

CALCULATED  EXAMPLES 

10 

1. 

The  NACA  64A006  With  C/4  Flap 

10 

1.1 

Steady  Flow  With  No  Flap  Deflection 

11 

1.2 

Steady  Flow  With  Deflected  Flap 

11 

1.3 

Unsteady  Flow 

16 

1.4 

Partial  Simulation  of  Viscous  Effects 

37 

1.5 

Simulation  of  Wind  Tunnel  Wall  Effects 

43 

2. 

Flow  Over  the  NACA  64A010  Airfoil 

58 

2.1 

Steady  Flow  at  Zero  Lift 

58 

2.2 

Steady  Flow  at  1. 0 Degree  Angle-of- Attack 

60 

2.3 

Indicial  Response  to  Step  Change  in  Angle-of- Attack 

60 

2.4 

Oscillatory  Unsteady  Flow 

66 

DISCUSSION 

81 

1. 

Correspondences  Between  Various  Calculations 

81 

2.  Correspondences  Between  Calculations  and  Experiments  85 
REFERENCES  87 


V 


LIST  OF  FIGURES 


Figure 

1 

2 

3 

4 

5 

6 

7 

8 

9 


Pressure  Distributions  on  NACA  64A006  at  Zero 
Angle-Of -Attack 

a.  Calculated 

b.  Experimental 

Pressure  Excursions  due  to  Flap  Deflection  in  Steady 
Flow.  NACA  64A006  at  Zero  Angle-Of-Attack 

a.  Mach  0.825 

b.  Mach  0.  85 

c.  Mach  0.  875 

Instantaneous  Pressure  Distributions  on  NACA  64A006 
with  Quarter-Chord  Oscillating  Flap,  Mach  0. 822,  Zero 
Angle-Of-Attack,  Flap  6 = 0 ±1.  5® , k = 0. 496 

Position  of  Upper  Surface  Shock,  NACA  64A006, 

Mach  0.822,  k = 0.496 

Instantaneous  Pressure  Distributions  on  NACA  64A006 
with  Quarter-Chord  Oscillating  Flap,  Mach  0. 854,  Zero 
Angle-Of-Attack,  Flap  6 = 0 ±1.0®,  k = 0.358 

Instantaneous  Pressure  Distributions  on  NACA  64A006 
with  Quarter-Chord  Oscillating  Flap,  Mach  0. 875, 

Zero  Angle-Of-Attack,  Flap  6 = 0 ±1.  0® , k = 0. 468 

Position  of  Upper  Surface  Shock,  NACA  64A006, 

Mach  0. 854,  k = 0. 358 

Position  of  Upper  Surface  Shock,  NACA  64A006, 

Mach  0. 875,  k = 0. 468 

Pressure  Excursions  Per  Radian  of  Flap  Deflection  in 
Unsteady  Flow.  NACA  64A006  at  Zero  Angle-Of-Attack, 
Mach  0.822,  k = 0.496 

a.  Amplitudes 

b.  Phase  Lead  Angles 


Page 


12 


13 


17 


20 


21 


24 

27 

27 


29 


Vi 


LIST  OF  FIGURES  (continued) 


Pressure  Excursions  Per  Radian  of  Flap  Deflection 
in  Unsteady  Flow.  NACA  64A006  at  Zero  Angle-Of- 
Attack,  Mach  0. 854,  k = 0.  358 

a.  Amplitudes 

b.  Phase  Lead  Angles 

Pressure  Excursions  Per  Radian  of  Flap  Deflection 
in  Unsteady  Flow.  NACA  64A006  at  Zero  Angle-Of 
Attack.  Mach  0.  875,  k = 0.468 

a.  Amplitudes 

b.  Phase  Lead  Angles 

Pressure  Distribution  on  NACA  64A006  at  Zero  Angle- 
Of-Attack,  Shock  Weakened  to  Simulate  Viscous 
Interaction,  Mach  0. 875 

Instantaneous  Pressure  Distributipns  on  NACA  64A006 
with  Quarter-Chord  Oscillating  Flap.  Shock  Weakened 
to  Simulate  Viscous  Interaction.  Mach  0. 875,  Zero 
Angle-Of-Attack,  Flap  6 = 0 ±1.  0° , k = 0.468 

Pressure  Distribution  on  NACA  64A006  at  Zero  Angle- 
Of-Attack  in  Ventilated  Wall  Tunnel,  Mach  0. 854 

Pressure  Distribution  on  NACA  64A006  at  Zero  Angle- 
Of-Attack  in  Ventilated  Wall  Tunnel,  Mach  0. 875 

Pressure  Excursions  Due  to  Flap  Deflection  in  Steady 
Flow.  NACA  64A006  at  Zero  Angle-Of-Attack  in  a 
Ventilated  Wall  Tunnel 

a.  Mach  0. 854 

b.  Mach  0. 875 

Instantaneous  Pressure  Distributions  on  NACA  64A006 
with  Quarter-Chord  Oscillating  Flap.  Zero  Angle-Of 
Attack  in  Ventilated  Wall  Tunnel.  Mach  0. 854,  Flap 
6 = 0 dhl.O',  k = 0.358 


LIST  OF  FIGURES  (Continued) 


Pressure  Excursions  per  Radian  of  Flap  Deflection  In 
Unsteady  Flow.  NACA  64A006  at  Zero  Angle -Of -Attack 
in  a Ventilated  Wall  Tunnel.  Mach  0.  854,  k = 0.358 

Pressure  Distributions  Along  a Line  1.  53  Chords  Above 
the  NACA  Airfoil.  Mach  0. 854,  Zero  Angle-Of-Attack, 
Quarter-Chord  Flap  Oscillating  at  k = 0. 358 

Instantaneous  Pressure  Distributions  on  NACA  64A006 
with  Quarter-Chord  Oscillating  Flap.  Zero  Angle-Of- 
Attack  in  Ventilated  Wall  Tunnel.  Mach  0. 875,  Flap 
a = 0 ±1.0”,  k = 0.468 

Pressure  Excursions  per  Radian  of  Flap  Deflection  in 
Unsteady  Flow.  NACA  64A006  at  Zero  Angle-Of-Attack 
in  a Ventilated  Wall  Tunnel.  Mach  0. 875,  k = 0.468 

a.  Amplitude 

b.  se  Lead  Angles 

istributions  on  NACA  64A010  at  Zero  Angle- 
ick  in  Mach  0. 80  Flow 

Pressure  Distributions  on  NACA  64A010  at  1.0  Degree 
Angle-Of-Attack  in  Mach  0. 80  Flow 

Calculated  Pressure  Distributions  on  NACA  64A010  in 
Mach  0.  80  Flow 

Pressure  Loading  on  NACA  64A010  Due  to  Step-Change 
in  Angle-Of-Attack  (a  = 0 ->1“ ) in  Mach  0. 80  Flow 

a.  Early  Stages 

b.  Later  Stages 

Transient  Lift  on  NACA  64A010  Due  to  Step-Change  in 
Angle-Of-Attack  (a  = 0 -*■1’ ) in  Mach  0. 80  Flow 

Transient  Pitching  Moment  on  NACA  64A010  Due  to 
Step-Change  in  Angle -Of- Attach  (a  = 0-^1”)  in  Mach  0.  80 
Flow.  Quarter-Chord  Reference  Axis,  Nose-Up  Positive 


r 


Figure 

28 


LIST  OF  FIGURES  (Concluded) 


Instantaneous  Pressure  Distributions  on  NACA  64A010  in 
Oscillatory  Plunging  Motion.  Mach  0.  80,  a = 0 ±1° , 
k = 0. 40 


Page 


67 


29 


30 


31 


32 


33 


Instantaneous  Pressure  Distributions  on  NACA  64A010  in 
Oscillatory  Pitching  Motion  About  a Quarter-Chord  Axis. 
Mach  0. 80,  a = 0 ±1” , k = 0. 40 

Instantaneous  Pressure  Distributions  on  NACA  64A010  in 
Oscillatory  Pitching  Motion  About  a Mid-Chord  Axis. 

Mach  0.  80,  0 = 0 ±1°  , k = 0.  50 

Position  of  Upper  Surface  Shock,  NACA  64A010  in  Mach 
0.  80  Unsteady  Flow 

Pressure  Excursions  due  to  1°  Angle-Of-Attack  Oscillations, 
NACA  64A010  in  Mach  0.  80  Flow 

a.  Amplitudes  in  Plunging  Motion 

b.  Phase  Lead  Angles  in  Plunging  Motion 

c.  Amplitudes  in  Pitching  Motion 

d.  Phase  Lead  Angles  in  Pitching  Motion 

Pressure  Distributions  for  NACA  64A006  at  Zero-Lift 
in  Mach  0.  85  Flow 


70 


73 


76 


77 


82 


Table 

1 


LIST  OF  TABLES 


Normalized  Forces  and  Moments  due  to  Flap  Deflection, 
NACA  64A006  with  Quarter-Chord  Flap,  First  Harmonic 
Values 

Force  and  Moments  in  1.  0 Degree  Oscillations  about  Zero 
Angle-Of-Attack,  NACA  64A010,  Mach  Number  0.80 


33 


79 


IX 


L 


SECTION  1 


INTRODUCTION 

Except  for  some  quite  restrictive  special  problems,  analysis  of  unsteady 
flows  over  airfoils  in  the  transonic  range  was  regarded  with  some  uneasiness 
for  a number  of  years.  Tolerance  for  the  use  of  numerical  solutions  on  such 
problems  has  grown,  and  many  efforts  to  provide  useful  transonic  computation 
procedures  have  been  carried  out.  Because  an  adequate  assessment  of  the 
probability  of  encountering  transonic  aeroelastic  instabilities  may  require 
checking  many  modes  over  a range  of  flight  conditions,  high  value  is  placed 
upon  keeping  the  aerodynamic  calculation  methods  simple.  Thus,  assessments 
of  the  perturbation  methods  by  comparison  of  solutions  done  by  such  methods 
with  calculations  done  by  methods  containing  fewer  restrictive  assumptions  are 
of  interest. 

The  work  being  reported  here  generated  numerical  solutions  to  a few 
problems  on  unsteady  transonic  flows  over  the  NACA  64A006  and  the 
NACA  64A010  airfoils.  These  examples  inelude  the  64A006  with  oscillating 
quarter-chord  flap  at  Mach  numbers  0. 822,  0. 854,  and  0.  875  and  the  64A010 
in  pitching  and  plunging  motions  at  Mach  number  0. 80.  Exploratory  calculations 
of  the  effect  on  the  oscillating  flap  results  of  wind  tunnel  walls  were  also  carried 
out. 

The  program  to  generate  these  solutions  is  based  upon  the  unsteady  Euler 
equations  in  conservation  form  and  the  finite- difference  procedure  uses  an 
explicit  shock-capturing  scheme.  The  choice  of  equations,  and  the  numerical 
method  preclude  the  method  being  applied  to  many  transonic  airfoil  problems 
and  the  results  cannot  show  flow  details  of  a scale  smaller  than  a few  mesh 
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widths.  Further,  one  should  keep  in  mind  that  these  are  exploratory  calcula- 
tions using  mesh  sizes,  means  for  assuring  stability  and  means  for  handling 
boundary  conditions  which  seemed  appropriate  for  the  ends  desired.  Changing 
the  computational  details  would  no  doubt  affect  the  results;  however,  little  has 
been  done  to  assess  the  effects  of  changing  the  details. 
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SECTION  II 


METHOD  OF  COMPUTATION 

The  unsteady  Euler  equations  in  conservation  form  are  solved  approxi- 
mately with  an  explicit  shocfc-capturing  scheme.  The  program  uses  features 

1,2 

quite  similar  to  those  described  in  earlier  woi’ks. 

1.  EQUATIONS 

The  computer  program  is  designed  to  obtain  numerical  approximations 
to  solutions  for  the  unsteady  Euler  equations  in  conservation  fomi: 


a(o)/3t  = -a(pu)/8x  - a(pv)/ay  (i) 

2 

a(pu)/at  = -a(pu  + p)/ax  - a(puv)/ay  (2) 

2 

a(pv)/8t  = -a(puv)/8x  - a(pv  +p)/8y  (3) 

8(E)/at  = - 8[u(E  + p)]  /8x  - a [V(E  + p)]  /8y  (4) 


Here,  the  usual  meanings  for  fluid  mechanics  problems  apply: 

p = fluid  density 
p = pressure 

u,v  = Cartesian  fluid  velocity  components 
x,y  = Cartesian  coordinates 
t = time 

y = adiabatic  index 
E = p/(Y  - 1)  + (P/2)(u^  + v^) 

(S 

= total  energy  per  unit  volume 


^Magnus,  R.  and  Yoshihara,  H. , "Inviscid  Transonic  Flow  over  Airfoils” 
AIAA  Journal,  Vol.  8,  No.  12,  pp.  2157-2162,  December  1970 
^Magnus,  R.  and  Yoshihara,  H.,  "Unsteady  Transonic  Flows  over  an  Airfoil" 
AIAA  Journal,  Vol.  13,  No.  12,  pp.  1622-1628,  December  1975 
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(6) 


The  coupled  system.  Equations  1 throu^  4,  has  the  form: 

8W/at  = -ar/ax  - ac/ay 

where  W,  F and  G are  four-element  vectors  whose  elements  may  be  discerned 
by  inspection. 

In  the  computer  program  the  system  represented  by  Equation  6 is  solved 
approximately  by  setting  up  arrays  of  about  5500  points  in  the  x-y  plane  at  which 
the  4 components  of  the  dependent  variable  W are  defined;  the  elements  of  F and 
G are  non-linear  functions  of  the  elements  of  W.  The  partial  derivatives  on  the 
right  hand  side  are  approximated  by  finite-differences  between  values  of  F and  G 
at  suitably  adjoining  mesh  nodes  near  the  point  being  worked  on  and  the  values  of 
W are  integrated,  approximately,  in  time  using  discrete  time  steps.  In  the  exist- 
ing computer  program,  airfoil  boundary  conditions  are  maintained  at  (on  the  order 
of)  100  of  the  mesh  nodes  distributed  along  a contour  which  would  coincide  with 
\ the  airfoil  if  the  airfoil  were  in  its  mean  position. 

I 2.  GRID  ARRANGEMENT 

I The  geometric  field  is  covered  with  a number  of  grids  having  a total  of 

j about  5500  nodes.  The  grids  of  the  outer  part  of  the  field  are  rectangular  with 

I uniform  spacing  in  the  chordwise  direction  and  non-untform  lateral  spacing 

I 

i v-stretched  spacing  outboard  of  | y|  = 0. 32  chord).  Within  a rectangle  circum- 

: cribed  about  the  airfoil  (-0. 26^  x s 1.34  chords  and  -0.32  « ys  0.32  chords) 
two  oTids  with  uniform  Ax  = 0. 04  chords  and  nonuniform  y - spacing  with  Ay 
^0. 04  chords  are  used.  These  two  grids  have  y - squashing  and  y - overlap  so 
that  there  are  nodes  of  these  systems  distributed  along  the  upper  and  lower  air- 
foil surfaces  except  on  the  blunt  part  of  the  nose.  The  two  systems  of  distorted 
lines  become  horizontal  and  coincident  a short  distance  behind  the  trailing  edge 
and  ahead  of  the  nose.  A portion  of  the  lenticular  gap  containing  and  conforming 
to  the  non-blunt  part  of  the  airfoil  is  covered  over  at  the  nose  by  a local  sheared 
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elliptic  system  conforming  to  the  blunt  part  of  the  airfoil  nose.  To  provide 
added  detail  of  shocks,  bands  of  Ax  = 0.  01  mesh  are  established  where  needed 
to  override  the  0. 04  chord  quasi-squares  immediately  surrounding  the  airfoil. 
The  shock-resolving  mesh  bands  are  0. 20  chord  wide  and  are  shifted  fore  and 
aft  in  0.  08  chord  steps  if  the  si  ock  moves  out  of  the  central  8 spaces  of  the 
band. 

The  central  block  of  rectangular  mesh  outboard  of  the  0. 04  chord  quasi- 
squares surrounding  the  airfoil  has  Ax  = 0. 08  chord  spacing  as  do  two  blocks 
of  mesh  ahead  of  and  behind  the  airfoil.  Further  i^stream  and  downstream 
there  are  blocks  with  Ax  = 0.16  and  Ax  = 0.32  chords;  the  field  is  thus 
extended  to  about  10  chords  ahead  of  and  behind  the  airfoil. 

The  spacing  in  each  coarser  rectangular  mesh  is  a power  of  two  larger 
than  the  spacing  in  each  finer  mesh,  and  lines  of  the  coarser  meshes  form  the 
perimeter  lines  and  some  of  the  interior  lines  of  the  finer  meshes.  This  con- 
formation of  grids  makes  the  data  exchanges  between  the  various  systems 
simple  and  also  simplifies  the  synchronization  of  the  calculations  that  pro- 
gress at  different  rates  in  the  various  grids. 

Waves  propagating  in  one  grid  may  be  partially  reflected  when  passing 
into  another  grid  of  different  coarseness  because  of  inevitable  differences  in 
the  dispersive  and  diffusive  truncation  errors  for  the  two  systems. 

3.  DIFFERENCE  SCHEME 

An  explicit  two-step,  finite-difference  scheme  of  the  Lax-Wendroff  type 

3 

was  used.  The  scheme  is  the  one  described  by  Thommen  and  has  been 
adapted  for  the  present  work  by  dropping  the  true  viscous  terms  and  substi- 
tuting a diffusion.  The  diffusive  damping  was  used  in  the  region  of  the  shock. 

3 

Thommen,  Hans  U. , "Numerical  Integration  of  the  Navier-Stokes  Equations" 
Journal  of  Applied  Mathematics  and  Physics,  (ZAMP)  Vol.  17,  Fasc.  3, 
pp.  369-384,  May  1966 
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In  calculating  steady  flows,  the  diffusive  damping  was  adjusted  according 
to  shock  strength  to  suppress  the  ragged  overshoots  at  shocks  which  are  a 
characteristic  of  Lax-Wendroff  difference  schemes.  On  unsteady  problems,  a 
single  (compromise)  damping  factor  was  usually  used  and  there  was  some 
tendency  for  weak  shocks  to  become  indistinct  because  of  unnecessarily  large 
damping. 

4.  SYNCHRONIZATION 

With  an  explicit  differencing  scheme,  the  calculation  is  stable  only  if  the 
time  step  is  limited.  Based  upon  typical  expected  flow  conditions  and  a know- 
ledge of  the  characteristics  of  the  scheme.  At  was  set  at  0. 35  Ax;  here,  a time 

unit  is  defined  as  C/a*  where  C is  the  airfoil  chord  and  a*  is  the  critical  speed 
00  00 

in  the  free  stream. 

In  advancing  the  solution  in  a finer  mesh  region  imbedded  in  a (say, 
twice  greater  mesh-size)  coarser  region,  two  time  steps  are  made  in  the 
finer  region  for  each  one  in  the  coarser.  Thereafter,  the  two  solutions  are 
adjusted  by  an  exchange  of  data.  The  solution  in  the  finer  mesh  is  assigned  to 
replace  the  coarser  mesh  solution  at  corresponding  points  in  the  interior  of  the 
finer  mesh,  and  the  solution  in  the  coarser  mesh  is  interpolated  onto  the  peri- 
meter of  the  finer  mesh  region.  Similar  principles  are  used  in  exchanging 
data  between  the  Cartesian  and  the  sheared  regions. 

5.  BOUNDARY  CONDITIONS 

The  problems  treated  here  all  deal  with  rigid  body  motion  of  the  profile 
or  flap.  The  mesh  system  is  independent  of  time  (except  for  the  fine-mesh 
patches  aroimd  shocks)  and  the  airfoil  boundary  conditions  are  satisfied  at 
those  mesh  nodes  which  lie  along  the  time-averaged  position  of  the  airfoil 
surface. 


When  tangency  boundary  conditions  are  to  be  satisfied  at  the  airfoil  sur- 
face the  advanced  values  of  the  dependent  variables  at  the  airfoil  nodes  are 
found  by  a sequential  process:  a)  A one-sided  difference  procedure  (obtaining 
data  only  from  the  exterior  of  the  airfoil)  is  used  to  obtain  tentative  values  for 
the  density,  energy  and  velocity  components,  b)  The  flow  velocity  component 
along  the  surface  normal  (whose  direction  may  be  a fxmction  of  time)  is 
adjusted  to  match  the  mechanical  velocity  of  the  surface  in  that  direction.  The 
adjustment  alters  the  normal  velocity,  density,  and  energy  in  the  manner  that 
they  would  be  altered  by  a plane  isentropic  wave  propagating  outward  from  the 
surface. 

Similar  wave  principles  are  used  to  match  pressures  and  flow  directions  ' 

along  the  upper  and  lower  sides  of  the  line  extending  aft  of  the  trailing  edge. 

The  geometric  trailing  edge  lies  between  mesh  nodes  and  no  further  specific  i 

steps  are  in  the  program  to  enforce  the  Kutta  condition. 

If  the  boundary  condition  to  be  enforced  were  to  be  the  assignment  of  a 

i 

particular  pressure,  the  wave  emitted  (substituting  for  step  (b)  in  the  j 

sequential  process  described  above)  would  be  chosen  to  cause  the  desired  |, 

pressure  and  the  "final"  velocity  normal  to  the  surface  would  be  a conse-  i 

quence  of  the  wave.  Suitable  modifications  to  this  wave  process  thus  allowed 
enforcement  of  free-jet  boimdary  conditions  and  idealized  conditions  at 
ventilated  wind  tunnel  walls. 

At  the  perimeter  of  the  computation  field  the  flow  variables  have  been 
held  fixed  at  some  pre-selected  values  during  the  solution  of  each  unsteady 
flow  problem.  When  steady  flow  problems  of  lifting  airfoils  were  solved,  the  ; 

circulation  on  the  perimeter  would  be  adjusted  periodically  to  match  the  lift 
(obtained  by  integrating  airfoil  pressures).  j 
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The  effects  of  some  of  the  approximations  in  satisfying  airfoil  and  field- 

4 

perimeter  boundary  conditions  were  studied  in  another  recent  work  and  found 
to  be  non-negligible  for  the  problems  studied.  The  mesh  system  used  in  the 
present  program  (with  y-stretching)  probably  lessens  the  errors  from  the 
inaccuracy  of  the  field  perimeter  boundary  conditions.  The  errors  due  to 
satisfying  unsteady  airfoil  boundary  conditions  at  fixed  nodes  rather  than  at 
nodes  on  the  translating  airfoil  surface  necessitates  a warning  to  the  reader 
that  the  "solutions"  generated  in  the  present  work  might  be  applicable  only  for 
vanishingly  small  oscillation  amplitudes. 

6.  COMPUTATIONAL  NOTES 

Geometric  information  on  each  problem  was  digested,  mesh  setup  opera- 
tions were  carried  out  and  interpolation  instructions  for  exchanges  of  data 
between  meshes  were  written  by  a preparatory  program.  The  time-integrations 
to  solve  the  unsteady  problems  were  carried  out  by  a main  program  which 
required  about  150K  (octal)  words  of  storage  on  a CDC  system.  The  programs 
are  coded  in  FORTRAN  extended  language  with  some  measures  to  take  advantage 
of  the  particular  capabilities  of  the  CDC  7600  central  processor.  An  FTN 
(OPT  = 2)  compilation  was  made  and,  generally,  the  problem  would  be  run  in 
stages  of  no  more  than  about  10  minutes  computing  for  each  stage.  This  would 
permit  inspection  of  progress  before  electing  to  continue.  A t)q)ical  problem, 
the  Mach  0.  875  flow  over  the  64A006  airfoil  with  flap  oscillating  at  reduced 
frequency  0. 468,  might  be  followed  for  3 1/2  oscillation  cycles  before  a 
(loosely  stated)  repeatable  oscillatory  response  was  obtained;  this  would 
require  about  2400  seconds  of  central  processor  time. 

4 

Magnus,  R.  J. , "Computational  Research  on  Inviscid,  Unsteady,  Transonic 
Flow  Over  Airfoils"  Convair  Rpt.  CASD/LVP  77-010,  January  1977 
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A few  exploratory  calculations  were  made  using  special  vector-algebra 
routines  (coded  for  the  CDC-7600  by  the  NASA-Ames  Computer  Fluid 
Dynamics  group)  in  the  principal  centered-finite-difference  scheme  of  the 
program.  There  was  little  increase  in  computational  speed,  probably 
because  the  computational  field  was  broken  into  small  blocks  of  mesh  having 
no  more  than  410  nodes  in  each  block  and,  in  fact,  the  vectors  operated  upon 
by  the  special  subroutines  n<  ^r  had  more  than  41  elements.  This  is  too 
small  a number  for  the  rouT;'  . to  be  efficient. 

Some  improvement  in  program  speed  was  accomplished  by  coarsening 
the  mesh  used  around  the  blunt  part  of  the  nose.  Because  the  mesh  around  the 
nose  is  smaller  than  anywhere  else  in  the  field,  the  explicit  difference  scheme 
must  make  many  time  steps  in  the  nose  region  for  each  time  step  in  the  0.  04 
chord  mesh  along  the  airfoil. 

The  program  had  been  constructed  specifically  to  calculate  flows  over 
the  64A006  with  oscillating  flap.  The  program  was  modified  repeatedly  to 
permit  calculations  over  pitching  and  plunging  airfoils,  assignment  of  pres- 
sures on  parts  of  the  airfoil  surface,  inclusion  of  wind  tunnel  walls,  etc. 

No  formal  documentation  on  the  program  has  been  prepared. 
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SECTION  III 


CALCULATED  EXAMPLES 


Calculations  were  made  on  two  airfoils  for  a small  number  of  operating 
conditions.  Because  of  the  expense  of  running  the  computer  program,  atten- 
tion has  been  concentrated  on  examples  which  would  display  transonic 
phenomena  of  interest  because  of  related  experimental  studies  or  to  furnish 
examples  for  comparison  with  calculations  done  by  other  methods. 

1.  THE  NACA  64A006  WITH  C/4  FLAP 

The  only  readily  available  modem  experimental  data  on  oscillatory 

5—1 0 

transonic  flow  over  airfoils  is  in  the  work  of  Tijdeman.  The  NACA  64A006 
airfoil  with  a quarter-chord  oscillating  flap  was  used  in  these  experiments  and 
this  configuration  was  the  subject  for  most  of  the  calculations  in  the  current 
research. 


Tijdeman,  H.  and  Bergh,  H. , "Analysis  of  Pressure  Distributions  Measured 
on  a Wing  With  Oscillating  Control  Surface  in  Two-Dimensional  High  Subsonic 
g and  Transonic  Flow"  Rpt.  NLR-TRF.253,  March  1967 
Tijdeman,  H.  and  Schippers,  P. , "Results  of  Pressure  Measurements  on  an 
Airfoil  with  Oscillating  Flap  in  Two-Dimensional  High  Subsonic  and  Tran- 
sonic Flow  (Zero  Incidence  and  Zero  Mean  Flap  Position),  Rpt.  NLR  TR 
^73078U,  July  1973 

Tijdeman,  H.  and  Schippers,  P. , "Results  of  Pressure  Measurements  on  a 
Lifting  Airfoil  with  Oscillating  Flap  in  Two-Dimensional  High  Subsonic  and 
gTransonic  Flow"  Rpt.  NLR  TR  73018L,  November  1974 
Tijdeman,  H. , "On  the  Motion  of  Shock  Waves  on  an  Airfoil  With  Oscillating 
gFlap  in  Two-Dimensional  Transonic  Flow"  Rpt.  NLR  TR  75038U,  March  1975 
Tijdeman,  H. , "High  Subsonic  and  Transonic  Effects  in  Unsteady  Aero- 
^gdynamlcs"  Rpt.  NLR  TR  75079U,  May  1975 
Tijdeman,  H. , "On  the  Unsteady  Aerodynamic  Characteristics  of  Oscillat- 
ing Airfoils  in  Two-Dimensional  Transonic  Flow"  Rpt.  NLR  MP  76003U, 
March  1976 
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1. 1 STEADY  FLOW  WITH  NO  FLAP  DEFLECTION 


Figure  la  shows  the  pressure  distribution  calculated  for  the  flow  over  a 

NACA  64A006  airfoil  at  zero  angle  of  attack  at  three  Mach  numbers.  Tijdeman's 

6,7 

experimental  pressure  distributions  are  reproduced  in  Figure  lb  for  com- 
parison. 

It  may  be  seen  that  there  are  appreciable  differences  between  the  calcu- 
lated and  experimental  results.  At  any  specified  Mach  number  the  calculated 
flows  are  more  supercritical  than  the  experimental  and  the  calculated  pressures 
on  the  aft  part  of  the  airfoil  are  more  positive  than  the  experimental  values. 
These  discrepancies  may  be  attributed  to  computational  inaccuracies  and  the 
influences  of  wind  tunnel  wall  interference  and  the  viscous  effects  which  are 
present  in  the  experiments  but  not  accounted  for  in  the  calculations. 

Inasmuch  as  the  magnitude  and  the  lag  of  unsteady  pressure  excursions 

(behind  the  flap  motion)  on  the  forward  parts  of  airfoils  with  oscillating  flaps 

depend  on  the  streamwise  lengths  and  the  heights  of  the  supersonic  regions 
5 8 9 10 

(see  Tijdeman  ’ * ' ) one  should  not  expect  to  find  good  agreement  between 

inviscid  unsteady  calculations  of  flow  over  an  airfoil  in  a free  stream  and  the 
unsteady  wind  timnel  experiments  at  any  specified  Mach  number  when  the  basic 
steady  flows  over  the  airfoil  differ  significantly. 

1.2  STEADY  FLOW  WITH  DEFLECTED  FLAP 

Pressure  distributions  calculated  at  Mach  numbers  0. 822,  0. 854,  and 

0. 875  for  the  64A006  with  quarter-chord  flap  deflected  by  1. 5,  1. 0 and  1.  0 

degrees  respectively  are  compared  to  the  experimental  results  of  Tijdeman  in 

Figures  2a  through  2c  which  present  the  normalized  pressure  excursions 

(calculated  as  (p,  - p )/6  flap)  at  a number  of  points  along  the  airfoil 

lower  upper 

chord.  In  each  case  the  calculated  pressure  excursions  exceed  the  measured 
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Figure  2a.  Pressure  Excursions  due  to  Flap  Deflection  in  Steady  Flow. 
NACA  64A006  at  Zero  Angle-Of-Attack  (continued) 


values  by  a large  factor  over  the  forward  portion  of  the  airfoil.  Traci  has 
also  calculated  this  flow  at  Mach  0.  85  and  0.  875,  and  except  for  some  parts  of 
the  airfoil  ahead  of  the  shock,  the  present  calculations  are  in  relatively  poor 
agreement  with  Traci's  results.  Normalizing  the  data  by  dividing  by  the 
magnitude  of  the  flap  deflection  unfortunately  does  not  lead  to  a uniform  pre- 
sentation of  the  loading  on  the  part  of  the  airfoil  traversed  by  the  displaced 
shock.  The  pressure  jump  across  the  shock  is  determined  by  the  upstream 

^^Traci,  R.  M. , Albano,  E.  D.  and  Farr,  J.  L.  Jr.,  "Small  Disturbance 
Transonic  Flows  About  Oscillating  Airfoils  and  Planar  Wings"  Rpt. 
AFFDL-TR-75-100,  August  1975 
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Pressure  Excursion,  (C 


Mach  0.85 
Figure  2b.  (Continued) 


flow  conditions  so  the  height  of  the  spike  in  the  vicinity  of  the  shock  tends 
toward  infinity  as  the  flap  deflection  is  reduced;  also  the  breadth  of  the  spike 
decreases. 


1.3  UNSTEADY  FLOW 

The  unsteady  flow  over  the  airfoil  was  calculated  at  the  following  three 
conditions: 


Mach  Number 


Reduced  Frequency 


Flap  Amplitude 


0.822 

0.854 


0.875 


0. 496 
0.358 
0.468 


Typical  pressure  distributions  at  six  phases  in  the  oscillation  cycle  for  these 
three  cases  have  been  presented  in  Figures  3,  5,  and  6.  The  three  cases  dis- 
play, respectively,  the  shock  motions  of  types  named  C,  B,  and  A by 

g 

Tijdeman;  see  Figures  4,  7,  and  8. 

The  three  imsteady  flow  problems  being  described  here  have  also  been 

12 

calculated  and  the  results  presented  by  Ballhaus  and  Goorjian.  Comparisons 

of  the  calculations  done  with  the  transonic  perturbation  program  and  those  done 

here  by  the  method  based  on  the  Euler  equations  were  presented  by  Ballhaus 
12 

and  Goorjian  for  the  cases  at  Mach  0. 854  and  0.  875.  Good  agreement 

between  the  results  using  the  two  methods  as  to  shock  location  was  demon- 

12 

strated  for  the  cases  at  Mach  0.  854  and  0. 875  by  Ballhaus  and  Goorjian. 


“^Ballhaus,  W.  F.  and  Goorjian,  P.  M. , "Implicit  Finite-Difference  Computa- 
tions of  Unsteady  Transonic  Flows  About  Airfoils,  Including  the  Treatment 
of  Irregular  Shock-Wave  Motions"  AIAA  Paper  77-205,  Los  Angeles, 
California,  January  1977 


16 


Figures  3a  & b.  Instantaneous  Pressure  Distributions  on  NACA  64A006  with  Quarter-Chord  Oscillating  Flap 
Mach  0.  822,  Zero  Angle-Of- Attack,  Flap  6 = 0 ±1. 5° , k = 0.  496  (Continued) 


0.4 


Figures  3c  & d.  (Continued) 


0.4 


0.7 


Figures  6e  & f.  Instantaneous  Pressure  Distributions  on  NACA  64A006  with  Quarter-Chord  Oscillating  Flap 
Mach  0.  875,  Zero  Angle-Of-Attack,  Flap  6 = 0 ±1.  0*,  k = 0.468  (Concluded) 


Fair  agreement  has  been  achieved  for  the  case  at  Mach  0. 822  but  the  current, 
Euler-based,  calculations  were  not  done  until  after  Reference  12  was  published. 

At  Mach  0. 822  the  flow  over  the  airfoil  with  undeflected  flap  is  sub- 
critical,  see  Figure  la.  With  the  flap  oscillating,  temporary  supersonic 
regions  with  weak  shocks  are  established  and  the  shocks,  in  sequence,  travel 
forward  to  the  nose  of  the  airfoil.  The  position  of  the  upper  surface  shock 
through  an  oscillation  cycle  is  shown  in  Figure  4,  compared  with  similar 
results  from  calculations  by  Ballhaus  and  Goorjlam^  Yu^^  has  also  calculated 
this  problem  but  used  a larger  flap  oscillation  amplitude. 

More  detail  on  the  propagation  of  the  shock  forward  on  the  airfoil  is 
furnished  in  the  pressure  distributions  shown  in  Figure  3g.  Limitations  on  the 
forward  movement  of  the  refined  "shock”  mesh  imposed  by  logic  in  the  current 
program  prevented  the  presentation  of  the  shock  profile  in  a 0. 01  chord  mesh 
if  the  shock  ran  forward  of  0. 20  chord.  Better  detail  of  the  shock  motion  on 
this  kind  of  problem  has  been  obtained  by  Yu  using  a transonic  perturbation 
program  with  shock  fitting. 

The  calculated  pressures  at  a number  of  loci  on  the  airfoil  have  been 

analyzed  to  determine  the  amplitude  and  phase  of  the  first  harmonic  responses 

to  the  sinusoidal  flap  motion.  These  data  are  presented  in  Figures  9,  10,  and 

6,7 

11  along  with  information  from  Tijdeman's  experiments. 

Comparing  the  amplitudes  of  the  calculated  pressure  excursions  at  the 

5,8,9,10 

three  Mach  numbers  it  is  evident,  as  has  been  observed  by  Tijdeman, 
that  at  the  higher  Mach  numbers  there  is  very  little  response  on  the  forward 

^®Yu,  N.  J.,  Seebass,  A.  R.,  and  Ballhaus,  W.  F. , "An  Implicit  Shock-Fitting 
Scheme  for  Unsteady  Transonic  Flow  Computations,"  AIAA  Paper  77-633, 
AIAA  3rd  Computational  Fluid  Dynamics  Conference,  Albuquerque, 

New  Mexico  June  27-28,  1977 
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Phase  Lead  Angles  x/C 

Figures  9a  & b.  Pressure  Excursions  Per  Radian  of  Flap  Deflection 
In  Unsteady  Flow.  NACA  64A006  at  Zero  Angle- 
Of-Attack.  Mach  0. 822,  k = 0. 496 
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Phase  Lead  Angles  x/C 

Figures  10a  & b.  Pressure  Excursions  Per  Radian  of  Flap  Deflection  in 

Unsteady  Flow.  NACA  64 A 006  at  Zero  Angle-Of-  Attack 
Mach  0.854,  k = 0.358 
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Figures  lla  & b.  Pressure  Excursions  Per  Radian  of  Flap  Deflection  in 
Unsteady  Flow.  NACA  64A006  at  Zero  Angle-Of-Attack 
Mach  0.  875,  k = 0.  468 


part  of  the  airfoil.  Signals  from  the  oscillating  flap  are  severly  delayed  and 
attenuated  in  getting  past  the  embedded  supersonic  regions. 

6 7 

Comparing  the  calculations  with  Tijdeman's  results  ’ and  looking  first 
at  the  flap  (X/C  > 0.75),  there  is,  on  the  whole,  fair  agreement  in  both  ampli- 
tude and  phase  of  the  pressure  excursions,  see  Figures  9,  10  and  11. 

On  the  part  of  the  airfoil  traversed  by  the  shocks,  the  calculated  excur- 
sions generally  exceed  the  measured  values  by  a factor  of  two  or  more.  This 
should  be  expected  because  the  calculation  should  show  a pressure  rise  at  the 
shock  approximating  "normal-shock"  recovery  whereas,  the  pressure 
immediately  downstream  of  the  shock  under  experimental  conditions  ought  to 
be  only  slightly  more  positive  than  the  pressure  for  Mach  one. 

On  the  part  of  the  airfoil  ahead  of  the  shock  (at  Mach  numbers  0.  854  and 
0. 875)  the  agreement  between  calculation  and  experiment  as  to  both  magnitude 
and  phase  of  the  pressure  excursions  is  rather  poor. 

16 

Ehlers  calculation  of  the  Mach  0. 854  case  and  the  present  results  are 
in  relatively  good  agreement  on  the  flap  but  differ  significantly  in  magnitude 
elsewhere.  The  phase  angles  agree  reasonably  well;  see  Figure  10. 

The  pressures  calculated  for  the  cases  at  the  three  Mach  numbers  men- 
tioned previously  were  integrated  to  find  the  airfoil  normal  force,  pitching 
moment  about  the  quarter  chord  (nose-up  positive)  and  the  flap  hinge  moment 
(also  considered  nose-up  positive).  These  have  been  listed  in  Table  1 for 


16 

Ehlers,  F.  Edward,  "A  Finite  Difference  Method  for  the  Solution  of  the 
Transonic  Flow  Around  Harmonically  Oscillating  Wings"  Rpt.  NASA 
CR-2257,  July  1974 
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TABLE  1.  NORMALIZED  FORCES  AND  MOMENTS  DUE  TO  FLAP  DEFLECTION, 

NACA  64A006  WITH  QUARTER-CHORD  FLAP,  FIRST  HARMONIC  VALUES 


J 


TABLE  1 (continued) 


Unstea^,  k = 0. 496 
Mach  Number  = 0. 822 


Experiment,  Tijdeman,  ’ 3.62  -24.  1.55  176.  0.0844  -167 

Mach  0. 858,  90  Hz 


comparison  with  Tijdeman's  experimental  data 

^ 11,12,13,16 

investigators. 


6,7 


and  calculations  by  other 


1.4  PARTIAL  SIMULATION  OF  VISCOUS  EFFECTS 

It  has  been  common  knowledge  since  the  earliest  experiments  on  transonic 
flow  over  airfoils  that  the  pressure  recovery  at  the  foot  of  a shock  at  the  down- 
stream border  of  an  imbedded  supersonic  zone  on  the  airfoil  is  less  than  normal- 
shock  recovery.  For  example,  T.  Theodorsen  (NACA  TN  1029,  March  1946) 
suggests  that  Mach  1.  0 is  attained  at  the  downstream  side  of  the  shock  and  refers 
to  earlier  works  by  H-S.  Tsien  and  A.  Fejer  (1944)  and  by  C.  L.  Dailey  (1943). 

In  the  absence  of  any  direct  evidence  on  the  pressure  just  aft  of  the  foot  of 

a shock  in  unsteady  flow,  it  was  assumed  that  the  pressure  would  be  close  to 

the  value  for  Mach  0.98.  As  a computational  device  for  weakening  the  shock  at 

the  airfoil  surface,  the  tangency  boundary  condition  was  relaxed  over  a limited 

portion  of  the  surface  (about  0. 1 chord)  downstream  of  the  moving  shock.  To 

establish  a procedure  for  the  shaping  and  placing  of  a lump-like  protuberance  to 

be  added  to  the  wall  aft  of  the  shock,  "reasonable"  pressure  profiles  based  upon 
6 7 

Tijdeman's  ’ measurements  were  prescribed  over  the  aft  part  of  the  64A006 
airfoil  and  the  computer  program  was  used  to  calculate  the  required  wall  shape 
modifications  for  shocks  of  several  strengths.  A generalized  description  of 
these  shapes  was  then  incorporated  in  the  unsteady  program. 
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Ballhaus,  W.  F.  and  Goorjian,  P.  M. , "Computation  of  Unsteady  Transonic 
Flows  by  the  Indicial  Method"  AlAA  Paper  77-447,  Dynamics  Specialist 
Conference,  San  Diego,  California,  March  24-25,  1977 
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In  calculating  an  oscillatory  flow  over  the  airfoil,  the  pressure  distribu- 
tion in  the  vicinity  of  the  shock  was  monitored  at  each  computational  pass  to 
establish  the  shock  location  and  strength  (based  principally  on  shock  upstream 
Mach  number  and  translation  speed  of  the  wave)  and  the  placement  and  shape  of 
the  lump  were  adjusted  to  cause  the  desired  less-than-normal  pressure  rise. 
The  computational  procedure  outlined  was  extremely  difficult  to  keep  stable. 

It  should  easily  be  recognized  that  the  lump  shape  would  have  little  con- 
nection with,  say,  the  displacement  thickness  profile  for  the  boundary  layer 
near  the  shock.  The  lump  shape  necessary  to  force  the  desired  pressure 
recovery  depends  upon  the  equations  used,  computational  scheme  and  mesh 
spacing.  This  can  be  (and  was)  easily  verified  by  prescribing  similar  airfoil 
pressures  and  calculating  the  lump  shapes  using  a program  based  upon  the 
transonic  perturbation  potential  with  non- conservative  differencing. 

Figure  12  shows  the  calculated  "viscous"  pressure  distribution  for  the 
64A006  at  zero  lift  in  a Mach  0.875  free  stream.  Representative  pressure 
distributions  at  6 phases  in  the  cycle  for  the  case  of  the  flap  oscillating  ±1. 0 
degree  at  reduced  frequency  k = 0.468  are  shown  in  Figure  13.  The  shock 
location  through  a cycle  is  shown  in  Figure  8;  compared  to  the  calculated 
inviscid  shock  motion,  the  shock  is  displaced  forward  and  the  response-wave 
form  is  substantially  altered.  The  deviation  of  the  wave-form  (of  shock 
location  through  the  oscillation  cycle)  from  the  shape  found  experimentally  by 
Tijdeman  most  likely  indicates  that  the  shock- weakening  computational  pro- 
cedure does  not  properly  simulate  the  consequences  of  the  real  (shock)- 
(boundary-layer)  interaction  on  the  pressure  field  aft  of  the  shock  and  on  the 
wave  propagation  along  the  airfoil  surface. 

The  oscillatory  force  and  moment  coefficients  for  this  case  with  partial 
simulation  of  viscous  effects  are  also  listed  in  Table  1 and  first  harmonies  of 
the  unsteady  pressures  are  shown  subsequently  in  Figure  21. 
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Figure  12.  Pressure  Distribution  on  NACA  64A006  at  Zero  Angle-Of-Attack, 
Shock  Weakened  to  Simulate  Viscous  Interaction,  Mach  0. 875 


The  programmed  procedures  for  simulating  viscous  effects  did  not  handle 
the  case  of  the  oscillatory  flow  over  the  64A006  in  a Mach  0.  854  stream  with 
flap  oscillating  ±1.  0 degree  at  k = 0. 358.  With  weaker  shocks,  the  lump  place- 
ment and  shaping  procedures  appeared  to  be  too  inexact  to  maintain  the  shock- 
lump  locked  on  the  downstream  end  of  the  siq)ersonic  region.  Revisions  to  the 
program  aimed  at  correcting  this  computational  instability  by  use  of  shock 
fitting  rather  than  shock  capturing  were  not  completed  satisfactorily. 
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Figures  13a  & b.  Instantaneous  Pressure  Distributions  on  NACA  64A006  with  Quarter-Chord  Oscillating  Flap 
Shock  Weakened  to  Simulate  Viscous  Interaction.  Mach  0.  875,  Zero  Angle-Of-Attack, 

Flap  6 =0*1.0°,  k 0. 468  (Continued) 


0.4 


wt  = 120' 


Figures  13e  & f.  Instantaneous  Pressure  Distributions  on  NACA  64A006  with  Quarter-Chord  Oscillating  Flap. 

Shock  Weakened  to  Simulate  Viscous  Interaction.  Mach  0.875,  Zero  Angle-Of -Attack, 

Flap  6 = 0 ±1.  0“  , k = 0.  468  (Concluded) 


1.  5 SIMULATION  OF  WIND  TUNNEL  WALL  EFFECTS 

The  calculations  described  in  Section  III  1. 1 through  1. 4 all  had  the 
perimeter  of  the  computation  field  at  least  10  chords  from  the  airfoil  to  try 
to  simulate  free-air  conditions.  In  this  section,  inviscid  calculations  of  the 
flow  over  the  64A006  in  a ventilated  wall  wind  tunnel  will  be  described. 

5 6 7 

Tijdeman's  experiments  on  the  unsteady  flow  over  the  64A006  ’ ’ 
airfoil  have  been  conducted  in  a slotted  wall  tuimel  with  a txmnel-height  of 
approximately  3.  06  chords.  The  slots  in  the  floor  and  ceiling  cover  about 
0. 10  of  those  areas.  The  walls  are  solid  ahead  of  the  test  section;  down- 
stream of  the  test  section  the  physical  arrangement  makes  it  appear  feasible 
to  apply  free- jet  boundary  conditions. 

Only  a limited  amount  of  data  on  the  behavior  of  the  pressures  near  the 

slotted  wall  in  unsteady  tests  in  this  particular  tunnel  has  been  made  avail- 
5 

able.  Therefore,  the  computational  model  of  the  slotted  part  of  the  wall  has 
been  based  upon  information  generated  in  calculations  performed  to  evaluate 
wall  effects  on  experiments  in  the  slotted  wall  tunnel  at  Wright- Patter  son 

17  17 

Air  Force  Base.  In  the  aforementioned  experiments  the  airfoil  pressures 
and  pressures  on  the  floor  and  ceiling  were  measured;  finite  difference  calcu- 
lations of  the  transonic  flowfields  using  these  data  determined  the  apparent 
vertical  velocity  at  the  tunnel  wall.  A rough  correlation  of  the  wall  velocity 
relations  in  these  experiments  has  been  made: 

u - 0.73V  - 0.17v  = 0 

X 

where  u and  v are  perturbations  from  free  stream  velocity  and  the  signs  are 
correct  for  the  lower  tunnel  wall. 

17 

Magnus,  R. , Yoshihara,  H,,  Lee,  D.,  and  Rogers,  L. , "Wall  Interference  in 
2D  Ventilated  Wind  Tunnels  at  High  Subsonic  Mach  Numbers"  Rpt.  AFFDL- 
TR-74-63,  June  1974 
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To  avoid  drastically  rearranging  the  computer  program,  the  calculations 
to  study  the  effect  of  the  wind  tunnel  walls  arbitrarily  have  used  solid  wall 
boundary  conditions  ahead  of  2. 6 chords  forward  of  the  model  and  free-jet 
conditions  aft  of  2.  6 chords  behind  the  model,  and  the  slotted  wall  con- 
dition along  the  intervening  part  of  walls. 

Calculations  on  the  NACA  64A006  with  quarter-chord  flap  in  a tunnel 
were  carried  out  at  Mach  numbers  0. 854  and  0. 875. 

The  zero-lift  steady  pressure  distributions  are  shown  in  Figures  14  and 
15  and  agree  much  closer  to  the  experimental  data  than  do  the  inviscid  free-air 
calculations,  refer  back  to  Figures  la  and  lb.  The  agreement  between  calcu- 
lation and  experiment  is  good  at  Mach  0. 875. 

In  a calculation  of  the  flow  over  the  64A006  at  Mach  0.  85  Traci^^  found 
that  including  a porous  wall  boundary  condition  made  the  flow  more  super- 
critical than  had  been  calculated  for  the  airfoil  in  a free  stream,  that  is,  the 
pressures  were  made  more  negative  on  the  midsection  of  the  airfoil  and  the 
shock  was  shifted  aft.  The  opposite  was  found  in  the  present  calculations,  no 
doubt,  because  of  a different  assumption  on  the  wall  characteristics. 

With  the  flap  deflected  1.  0 degree  in  steady  flow,  the  pressure  excur- 
sions on  the  forward  part  of  the  airfoil  are  greatly  reduced  from  the  values 
noted  for  the  calculations  in  a free-stream;  see  Figures  16  and  2.  In  this 
respect,  the  present  calculations  agree  with  the  trend  found  by  Tracl^^  whose 
results  are  also  shown  in  Figure  16a. 

In  Mach  0. 854  flow.  Figures  17  and  18,  including  wind  tunnel  wall  effect 

in  the  calculations  of  the  flow  with  the  flap  oscillating  produces  relatively  good 

6,7 

agreement  with  Tijdeman's  experimental  data  in  most  respects.  The 
calculation  still  predicts  larger  pressure  excursions  on  the  part  of  the  airfoil 
traversed  by  the  moving  shock  than  are  observed  in  the  experiments;  this  is 
inevitable  unless  the  shock  is  deliberately  weakened. 
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Figure  14.  Pressure  Distribution  on  NACA  64A006  Figure  15.  Pressure  Distribution  on  NACA  64A006  at 

at  Zero  Angle-Of-Attack  in  Ventilated  Zero  Angle-Of-Attack  in  Ventilated  Wall 

Wall  Tunnel,  Mach  0.854  Tunnel,  Mach  0.  875 


Figures  16a  & b.  Pressure  Excursions  Due  to  Flap  Deflection  in  Steady  Flow.  NACA  64A006 
at  Zero  Angle-Of-Attack  in  a Ventilated  Wall  Tunnel 
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Figures  17a  & b.  Instantaneous  Pressure  Distributions  on  NACA64A006  with  Quarter-Chord  Oscillating 
Flap.  Zero  Angle-Of-Attack  in  Ventilated  Wall  Tunnel.  Mach  0.854,  Flap  6 = 0 ±1.0' 
k = 0.  358  (Continued) 
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Detail  of  Lower  Surface  Pressures, 
ut  = 150*  to  180*  x/C 

Figure  17g.  Instantaneous  Pressure  Distributions  on  NACA  64A006  with  Quarter-Chord 
Oscillating  Flap.  Zero  Angle-Of-Attack  in  Ventilated  Wall  Tunnel.  Mach 
0.  854,  Flap  6 = 0 ±1.  0°  , k = 0.  358  (Concluded) 

The  distributions  of  pressure  along  a line  1.  53  chords  above  the  airfoil 
are  shown  in  Figure  19  for  calculations  in  a ventilated  wall  tunnel  and  in  a free 
stream.  The  wall  shifts  the  peak  negative  pressure  aft  and  reduces  the  magni- 
tude of  the  peak.  The  amplitude  of  the  unsteady  part  of  the  pressures  has  a 
tendency  to  be  large  when  the  mean  pressure  is  large;  compare  Figures  l9a 
and  b.  The  phase  lead  angles  are  affected  somewhat  also. 

In  Mach  0. 875  flow  with  the  flap  oscillating  sinusoidally  the  pressure 
excursions  over  the  forward  part  of  the  airfoil  are  greatly  reduced  in  magni- 
tude from  their  values  in  free  air,  see  Figures  20,  21,  and  11.  In  fact,  the 

pressure  excursions  calculated  for  the  forward  part  of  the  airfoil  are  lower  in 

6 7 

magnitude  than  the  experimental  values  of  Tijdeman.  ’ The  calculated  phase 
angles  for  the  pressure  excursions  are  brought  markedly  closer  to  the 
experimental  values  by  Including  wind  tunnel  wall  effect.  These  differences 
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Figures  18a  & b.  Pressure  Excursions  per  Radian  of  Flap  Deflection  in 
Unsteady  Flow.  NACA  64A006  at  Zero  Angle-Of-Attack 
In  a Ventilated  Wall  Tunnel.  Mach  0. 854,  k = 0. 358 


Mean  Pressures 


Figure  19a.  Pressure  Distributions  Along  a Line  1.53  Chords  Above  the  NACA 
Airfoil.  Mach  0.854,  Zero  Angle-Of-Attack,  Quarter-Chord 
Flap  Oscillating  at  k = 0. 358  (Continued) 

on  the  forward  part  of  the  airfoil  have  only  a small  effect  on  the  airfoil  oscilla- 
tory forces  and  moments  at  Mach  0.  875  since  the  excursions  of  pressure  on 
the  forward  part  of  the  airfoil  are  negligible  compared  to  those  on  the  flap  or  in 
the  region  traversed  by  the  moving  shock. 

The  positions  of  the  shocks  in  unsteady  flow  in  the  ventilated-wall  tunnel 
are  shown  in  Figures  7 and  8. 

g 

At  Mach  0.  854  in  the  tunnel,  the  shock  motion  seems  to  be  of  Tijdeman's 
Type  C.  This  is  evident  in  the  pressure  distributions  of  Figure  17  but  the  pro- 
gram logic  did  not  insert  mesh  to  provide  detail  of  the  pressure  rise  in  this 
shock  as  it  ran  forward  in  subsonic  flow  toward  the  nose.  In  free  air  at 
Mach  0. 854  the  shock  motion  seems  to  be  of  Type  6. 
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Figures  19b  & c. 
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Pressure  Distributions  Along  a Line  1. 53  Chords  Above  the 
NACA  Airfoil.  Mach  0. 854,  Zero  Angle-Of-Attack, 

Quarter  Chord  Flap  Oscillating  at  k « 0. 358  (Concluded) 
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Figures  20a  & b.  Instantaneous  Pressure  Distributions  on  NACA  64A006  with  Quarter-Chord  Oscillating  Flap. 

Zero  Angle-Of-Attack  in  Ventilated  Wall  Tunnel.  Mach  0, 875,  Flap  6 = 0 ±1. 0* , k = 0. 468  (Cont) 


35.0 
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Figure  21a.  Pressure  Excursions  per  Radian  of  Flap  Deflection  in  Unsteady 
Flow.  NACA  64A006  at  Zero  Angle-Of-Attack  in  a Ventilated 
Wall  Tunnel.  Mach  0. 875,  k = 0. 468  (Continued) 

14 

Ballhaus  has  found  by  comparing  calculations  in  free  air  and  in  a 
restricted  stream  with  free- jet  boundtmies,  that  shock  motion  pattern  could 
be  altered  from  Type  A in  free  air  to  Type  B in  the  free-jet. 

14 

Ballhaus,  W.  G.  and  Goorjian,  P.  M. , "Efficient  Solution  of  Unsteady  Trans- 
sonic Flows  About  Airfoils"  44th  AGARD  Structures  and  Materials  Panel 
Meeting,  Lisbon,  18  April  1977 
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Figure  21b.  Pressure  Excursions  per  Radian  of  Flap  Deflection  in  Unsteady 
Flow.  NACA  64A006  at  Zero  Angle-Of-Attack  in  a Ventilated 
Wall  Tunnel.  Mach  0.875,  k = 0.468  (Concluded) 


2.  FLOW  OVER  THE  NACA  64A010  AIRFOIL 

Flowfields  around  the  NACA  64A010  airfoil  at  Mach  0. 80  have  been  calcu- 
lated for  four  unsteady  cases.  This  airfoil  is  to  be  the  initial  model  in  a series 
of  unsteady  transonic  experiments  to  be  carried  out  at  NASA  Ames  Research 
Center. 


2. 1 STEADY  FLOW  AT  ZERO  UFT 

The  pressure  distribution  calculated  at  zero  angle-of-attack  for  the 
symmetric  airfoil  is  shown  in  Figure  22.  Similar  calculations  by  Ballhaus 

and  Goorjian  using  a program  based  on  the  transonic  peiturbation  potential 

12  18 
equation  and  by  Steger  using  a program  based  on  the  Euler  equations  are 

. 6 

shown.  Experimental  data  from  Stivers  (unpublished,  Re  = 4 x 10  , natural 
transition)  are  also  presented. 

^^Steger,  J.  L. , "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  With  Application  to  Airfoils"  AIAA  Paper  77-665, 
AIAA  10th  Fluid  and  Plasma  Dynamics  Conference,  Albuquerque, 

New  Mexico,  June  27-29,  1977 
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2.  2 STEADY  FLOW  AT  1. 0 DEGREE  ANGLE -OF -ATTACK 
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At  1. 0 degree  angle -of- attack  in  Mach  0. 80  flow  the  calculated  pressure 
distribution  shows  very  little  lift  being  generated  on  the  aft  part  of  the  airfoil. 
Figure  23.  The  present  calculations  agree  well  with  those  of  Steger  which  were 
made  with  a program  (similarly)  based  upon  conservation  equations  but  show 
less  lift  than  that  calculated  by  Ballhaus  and  Goorjian.  To  assist  in  interpret- 
ing the  unsteady  pressure  distributions,  the  steady  flow  calculations  are  also 
presented  as  p/Pq  in  Figure  24. 

The  pressure  excursions  due  to  1.0  degree  angle-of- attack,  calculated 

as  (p,  - p )/2a!  are  shown  subsequently  in  Figure  32. 

lower  upper 

2. 3 INDICIAL  RESPONSE  TO  STEP  CHANGE  IN  ANGLE-OF -ATTACK 

Starting  from  a flow  solution  for  zero  lift,  an  increment  in  cross -flow 
velocity  equal  to  0. 01745  was  added  to  the  velocity  at  each  mesh  node  at 
zero  time. 


x/C 

Figure  24.  Calculated  Pressure  Distributions  on  NACA  64A010  in 
Mach  0.  80  Flow. 
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Early  Stages 

Figure  25a.  Pressure  Loading  on  NACA  64A010  Due  to  Step-Change  in  Angle- 
Of -Attack  (Q!  = 0 — >1* ) in  Mach  0.  80  Flow  (Continued) 

Chordwise  distributions  of  the  increments  in  pressure  loading  on  the  air- 
foil at  several  instants  in  the  first  chord  and  a half  of  travel  are  shown  in 
Figure  25.  During  this  time  period,  the  rectangular  initial  loading  changes 
toward  a loading  pattern  peaked  at  the  nose  and  the  loading  due  to  shock  move- 
ment becomes  apparent.  Here,  a time  unit  is  C/U^ . 

The  unsteady  lift  and  pitching  moment  about  the  quarter  chord  are  shown 
as  functions  of  time.  Figures  26  and  27,  and  are  compared  with  the  results  of 
similar  calculations  by  Ballhaus  and  Goorjian  and  by  Steger. 

The  initial  parts  of  the  transients  calculated  here  differ  greatly  from 
those  calculated  by  Ballhaus  and  Goorjian  using  the  perturbation  equations. 
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Figures  26a  & b.  Transient  Lift  on  NACA  64A010  Due  to  Step-Change  in 
Angle-Of-Attack  (or  = 0 -^1®)  in  Mach  0.80  Flow 
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The  sudden  imposition  of  the  cross-flow  induces  a relatively  uniform 
jump  in  pressure  on  the  lower  airfoil  surface  and  a similar  fall  in  pressure  on 
the  upper  airfoil  surface;  see  curves  labelled  t = . 0115  in  Figure  25a.  The 
lower  surface  compression  and  the  upper  surface  rarefaction  propagate  their 
influence  to  the  respective  opposite  surfaces  around  both  the  leading  and  trailing 
edges. 

There  is  appreciable  attenuation  of  each  of  the  primary  disturbances  in 
spreading  around  the  leading  edge  to  the  opposite  surface.  The  partial  can- 
cellation of  the  primary  waves  around  the  leading  edge  is  swept  aft  at  speed 
(u  + a)  causing  some  drop  in  lift.  These  adjustments  sweep  to  the  trailing  edge 
in  about  0.4  time  units.  Likewise,  the  mutual  cancellation  of  the  primary  dis- 
turbances around  the  trailing  edge,  which  spreads  forward  at  speed  (a  - u),  is 
well  Illustrated  in  Figure  25.  During  the  entire  period  covered  in  Figure  25 
the  loading  at  the  leading  edge  builds  up  and  spreads  aft. 

The  wave  systems  mentioned  above  are  all  features  present  in  the 

19 

linearized  analysis  of  Heaslet  and  Lomax.  Because  the  airfoil  is  not  a flat 
plate  starting  with  uniform  flow  as  assumed  in  the  linear  analysis,  the  various 
waves  of  adjustment  interact  with  the  non-uniform  flow  around  the  finite-thick- 
ness aiiioil.  Thus,  loading  spikes  develop  due  to  movement  of  the  shocks 
(Figure  25b)  and  waves  travelling  forward  from  the  trailing  edge  are  severly 
attenuated  and  delayed  by  having  to  detour  outboard  of  the  shocks  closing  the 
imbedded  supersonic  regions. 

Both  the  lift  and  pitching  moment  decrease  in  magnitude,  Figures  26a 
and  27a,  to  minimums  at  about  0. 4 time  units  after  their  initial  jumps,  about 


Heaslet,  Max.  A.  and  Lomax,  Harvard,  "The  Application  of  Green's  Theorem 
to  the  Solution  of  Boundary-Value  Problems  in  Linearized  Si;g)ersonic  Wing 
Theory"  National  Advisory  Committee  for  Aeronautics  Technical  Note  No. 
1767,  April  1949 
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the  time  for  information  to  propagate  aft  along  the  airfoil  one  chord  length. 
The  pitching  moment  holds  steady  at  about  2/3  its  "final"  value  for  the  period 
between  (roughly)  4 and  7 time  units  after  the  step  change  in  incidence;  a 
similar  pause  is  present  in  the  moment  buildup  calculated  by  Ballhaus  and 
Goorjian  and  by  Steger,  Figure  27a. 

2.4  OSCILLATORY  UNSTEADY  FLOW 

The  oscillatory  unsteady  flow  over  the  airfoil  at  Mach  0. 80  was  calcu- 
lated for  three  oscillation  modes: 


Mode 

Axis 

Incidence 

Reduced 

F reauency 

Pitching 

0. 50  chord 

0 dd.0“ 

0.  50 

Pitching 

0. 25  chord 

0 ±1.0" 

0.40 

Plunging 

0 ±1.  0*  (Induced  by 
Motion) 

0. 40 

Instantaneous  pressure  distributions  at  six  positions  in  the  oscillation  cycle 
for  the  plunging  motion  are  shown  in  Figure  28.  The  pressure  distributions 
for  the  cases  with  pitching  motion  are  similar  in  character  but,  of  course, 
differ  somewhat  in  their  detailed  values.  Figures  29,  30.  Figure  31  shows 
the  upper  surface  shock  location  through  an  oscillation  for  the  three  cases 
calculated. 

The  amplitudes  and  phase  angles  for  the  pressure  excursions  for  the 
three  cases  are  shown  in  Figure  32. 

The  oscillatory  force  and  pitching  moment  about  the  quarter  chord  (nose- 
up  is  positive)  are  listed  in  Table  2.  In  steady  flow  at  1. 0 degree  angle-of- 
attack  the  present  results  agree  satisfactorily  with  those  of  Steger.  The  lift 
calculated  by  Ballhaus  and  Goorjian  is  about  1.3  times  the  average  of  the  other 
two  results.  Similarly,  the  pitching  moment  calculated  by  the  perturbation 
method  is  1. 6 times  the  average  of  the  other  two  results. 


Figures  28a  & b.  Instantaneous  Pressure  Distributions  on  NACA  64A010  in  Oscillatory  Plunging 
Motion.  Mach  0.  80,  a = 0 ±1*  , k = 0.  40  (Continued) 


Figures  29a  & b.  Instantaneous  Pressure  Distributions  on  NACA  64A010  in  Oscillatory  Pitching 

Motion  About  a Quarter-Chord  Axis.  Mach  0.80,  Oil",  k = 0.40  (Continued) 


Figures  29c  & d.  (Continued) 


Figures  29e  & f.  Instantaneous  Pressure  Distributions  on  NACA  64A010  in  Oscillatory  Pitching 

Motion  About  a Quarter-Chord  Axis.  Mach  0.  80,  a = 0±1°  , k = 0.  40  (Concluded) 


Figures  30a  & b.  Instantaneous  Pressure  Distributions  on  NACA  64A010  in  Oscillatory  Pitching  Motion 
About  a Mid-Chord  Axis.  Mach  0.80,  a - 0 ±1“,  k = 0. 50  (Continued) 


About  a Mid-Chord  Axis.  Mach  0.80,  o = 0 ±1* , k = 0.  50  (Concluded) 


wt 

Figure  31.  Position  of  Upper  Surface  Shock,  NACA  64A010  in  Mach  0.80 
Unsteady  Flow 

In  plunging  unsteady  motion,  the  magnitude  of  the  lift  calculated  by 

Ballhaus  and  Goorjian  by  the  indicial  method  using  the  transonic  perturbation 

equations  is  also  higher  than  the  other  calculations.  The  phase  angles  for  the 

lift  show  a spread  of  11  degrees  for  the  various  calculations.  The  unsteady 

pitching  moments  for  the  plunging  problem  vary  in  magnitude  from  .0105  to 

. 0186.  Among  the  various  calculations  there  is  an  overall  spread  in  the  phase 

angles  for  the  pitching  moment  of  about  23  degrees.  The  calculations  having 

the  greatest  difference  in  magnitude  for  the  xmsteady  pitching  moment  are  the 

present  and  an  analogous  calculation  with  airfoil  boundary  conditions  satisfied 

4 

at  nodes  on  the  plunging  airfoil  rather  than  at  a time-averaged  position  of  the 
airfoil. 

The  differences  between  the  results  of  the  two  pitching  problems  calcu- 
lated are  relatively  minor  (Table  2). 
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Table  2.  FORCES  AND  MOMENTS  IN  1. 0 DEGREE  OSCILLATIONS  ABOUT  ZERO 
ANGLE-OF -ATTACK,  NACA  64A010,  MACH  NUMBER  0.  80 
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SECTION  IV 
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DISCUSSION 

The  calculations  completed  ought  to  have  two  uses;  first,  as  examples  of 
what  might  be  expected  from  a straightforward  munerical  solution  of  the 
inviscid  unsteady  Euler  equations  as  compared  to  other  solutions  and,  second, 
as  an  aid  in  interpreting  what  has  been  found  in  experiments. 

1.  CORRESPONDENCES  BETWEEN  VARIOUS  CALCULATIONS 

When  checking  for  agreement  between  results  of  various  unsteady  calcu- 
lations one  should  expect  differences  to  occur  because  of  the  diversity  in 
formulations  of  the  problem  and  in  the  mechanizations  for  approximate 
numerical  solution. 

Most  investigators  have  calculated  steady  flows  around  the  airfoil  being 
studied  to  establish  a limit  of  zero  reduced  frequency  or  as  a basic  step  In  the 
solution  sequence.  It  can  be  seen.  Figure  33,  that  there  is  considerable 
spread  in  the  results  of  solving  what  might  be  expected  to  be  a simple  flow 
problem,  namely  flow  past  a thin  airfoil  at  zero  lift.  There  is  only  fair  agree- 
ment on  obvious  features  which  ought  to  be  important,  namely  the  size  and 
placement  of  the  supersonic  region. 

When  the  unsteady  problem  imder  stu(fy  is  the  oscillating  flap,  widely 
different  results  could  be  expected  in  the  responses  if  the  basic  flows  had 
different-sized  supersonic  regions  because  of  the  different  delays  and  attenua- 
tions of  the  signals  from  the  moving  flap  reaching  the  front  of  the  airfoil  via 
paths  which  (certainly)  extend  outboard  of  the  imbedded  supersonic  regions. 
When  the  problem  being  studied  is  pitching  or  plunging,  the  interplay  of  signals 
originating  from  changes  and  attitude  or  velocity  of  the  elements  on  the  airfoil 
surface  as  a whole  is  more  complicated  than  for  the  oscillating  flap  problem. 
Still,  the  delays  in  propagation  of  information  forward  from  the  trailing  edge 
are  important  in  determining  the  load  on  the  complete  airfoil. 
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Figure  33.  Pressure  Distributions  for  NACA  64A006  at  Zero-Lift  in  Mach  0.  85 
Flow 

Even  if  there  is  little  prospect  of  absolute  agreement  between  methods  it 
would  be  of  value  to  determine  whether  the  methods  predict  the  same  trends  as 
parameters  (such  as  Mach  number,  reduced  frequency  and  amplitude  of  oscilla- 
tion) are  varied.  There  might  be  agreements  with  some  bias-like  shifts.  That 
is,  increasing  the  Mach  number  slightly  for  method  "A"  might  make  it  agree 
better  with  method  "B"  because  the  two  then  would  show  flows  with  roughly 
the  same  basic  supercriticality.  Likewise,  method  "C"  might  consistently 
show  less  pitching  moment  than  method  "D"  because  in  "C"  a coarser  mesh 
was  employed  and  the  shock  capturing  scheme  places  the  shock  consistently 
further  forward  on  the  airfoil  than  for  method  "D". 

Unfortunately,  of  the  literature  examined,  there  have  been  only  two 

12  13  11 

research  teams,  namely  Ballhaus  and  Goorjian  * and  Traci,  et  al,  who 
20 

Isogai,  Koji,  "Calculation  of  Unsteady  Transonic  Flow  Over  Oscillating 
Airfoils  Using  the  Full  Potential  Equation"  AIAA  Paper  77-448, 

AlAA  Dynamics  Specialist  Conference,  San  Diego,  California 
March  24-25,  1977 
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have  calculated  systematic  sequences  of  cases.  The  number  of  cases  calcu- 
lated in  the  present  research  is  too  small  and  the  parameters  have  not  been 
varied  systematically  enough  to  demonstrate  an  orderly  dependence  of  response 
characteristics  on  the  problem  parameters. 

Where  the  parameters  for  the  calculations  were  deliberately  picked  so 
that  direct  comparisons  could  be  made  between  the  present  method  and  that  of 
Ballhaus  and  Goorjian,  ’ the  relative  values  predicted  by  the  two  methods 
have  been  puzzling.  Referring  to  Table  1,  at  Mach  0. 822  the  present  method 
predicts  a lower  oscillating  lift  than  the  method  of  Ballhaus  and  Goorjian;  at 
Mach  0. 875,  the  opposite  is  true.  Also,  the  disagreement  between  the  two 
styles  of  calculation  seems  greatest  at  the  lowest  Mach  number  examined. 

On  the  problem  of  Mach  0. 80  flow  over  the  64A010  in  plunging  oscillations 

or  in  a step  change  in  angle-of-attack,  the  disagreement  between  the  present 

13 

results  and  those  of  Ballhaus  and  Goorjian  seems  to  stem  from  the  non- 
agreement as  to  the  steady  flows. 

Where  a direct  comparison  has  been  made  between  the  present  method 

(modified  to  satisfy  boundary  conditions  at  nodes  which  move  with  the  plunging 
4 18 

airfoil)  and  that  of  Steger,  see  Table  2,  the  agreement  has  been  relatively 

good.  The  two  methods  are  similar  in  solving  a coupled  system  of  four 

vmsteady  conservation  equations  with  shock  capturing  schemes.  The  programs 
18 

differ  in  that  Steger  uses  a single  automatically  generated  mesh  and  an 
implicit  scheme  and  is  coded  to  solve  viscous  flows  usii^  the  Navier  Stokes 
equations  with  a turbulence  model;  the  program  used  here  has  used  patches  of 
different  sized  mesh  to  provide  the  resolution,  an  explicit  scheme,  and  is  coded 


for  the  inviscid  Euler  equations.  The  effects  of  two  approximations  made  in 
the  immediate  antecedent  to  the  present  program,  handling  of  boundary  con- 


ditions at  the  airfoil  surface  and  the  assignment  of  flow  properties  at  the  peri- 

4 

meter  of  the  computation  field  were  studied  recently.  The  study  demonstrated 
that  modest  changes  in  the  magnitudes  and  phase  angles  for  lift  and  pitching 
moment  response  could  be  expected  if  the  calculation  were  to  be  performed  in 
a truly  unbounded  computation  field  rather  than  the  9. 6 chord  square  which  had 
been  used  in  that  study.  The  present  calculations  have  used  a stretched  y-mesh 
in  an  attempt  to  lessen  the  effect  of  Imperfection  in  applying  the  outer  boundary 
condition  but  the  errors  occasioned  by  this  practice  have  not  been  evaluated. 
Applying  the  airfoil  boundary  conditions  at  fixed  nodes  located  at  the  time- 
averaged  position  of  the  oscillating  airfoil  (which  is  what  was  done  in  all  calcu- 
lations in  the  present  study)  increased  the  magnitude  and  increased  the  lag 
angle  for  the  pitching  moment  by  about  16  degrees  on  the  one  example  studied 
see  Table  2,  Section  III  2.  of  the  present  report.  Ostensibly  these  particular 
errors  should  be  proportional  to  the  linear  velocities  at  airfoil  surface 
elements  and  are,  therefore,  dependent  on  oscillation  amplitude.  Thus,  it 
would  be  best  to  regard  the  results  of  the  calculations  in  the  present  report  as 
being  applicable  to  problems  with  vanishingly  small  amplitude.  The  only  pro- 
gram, of  those  examined  here,  incorporating  a combination  of  equations  and 

boundary  conditions  which  seems  suitable  for  problems  involving  thick  airfoils 

18 

in  finite  amplitude  oscillations  is  that  of  Steger. 

In  seeking  data  for  comparisons  of  the  present  results  with  those  of  other 
investigators  it  was  apparent  that  there  is  not  much  standardization  on  presen- 
tation of  data.  Thus,  we  can  find  reduced  frequency  defined  more  than  one 
way,  oscillatory  data  presented  as  real  and  imaginary  components  and  as 
magnitudes  and  phase  angles,  diverse  non-dimensionallzations  to  obtain 
coefficients,  pitching  moments  about  the  nose,  the  quarter-chord  and  mid- 
chord. This  is  not  a suggestion  that  all  investigators  should  solve  the  same 


problems  and  present  the  results  in  the  same  fashion.  Certainly  programs 
which  are  based  on  particular  abridged  equations  will  be  useful  on  problems 
which  do  not  violate  the  basic  assumptions  and  new  programs  will  be  needed  to 
work  on  problems  which,  because  of  their  peculiar  features  or  severity,  cannot 
be  handled  by  present  methods. 

2.  CORRESPONDENCES  BETWEEN  CALCULATIONS  AND  EXPERIMENTS 

The  ultimate  purpose  of  calculations  should  be  the  prediction  of  the 
behavior  of  oscillating  wings  at  flight  Reynolds  numbers.  Considering  cost  and 
other  feasibility  factors,  it  is  unlikely  that  experimental  data  under  fli^t  con- 
ditions will  be  available  soon.  Thus,  the  programs  ought  to  be  capable  of 
duplicating  findings  in  wind  tunnel  tests  and  then  extendable  by  logical  changes 
to  predict  behavior  under  flight  conditions.  The  changes  should  account  for 
Reynolds  number  differences  between  flight  and  wind  tunnel  conditions  and  for 
variations  in  external  influences,  that  is  three-dimensional  geometry  in  free- 
air  rather  than  two-dimensional  or  other  idealized  geometiy  in  a channel  with 
ventilated  walls. 

21 

According  to  Kacprzynski,  the  combined  uncertainties  as  to  the  precise 
effects  of  wind-tunnel  wall  interference  and  the  effects  of  viscosity  may  result 
in  20%  error  between  theory  and  experiment  on  studies  of  steady  transonic 
flows  over  modem  supercritical  airfoils.  It  is  probably  unreasonable  to 
believe  that  much  better  agreement  will  be  achieved  in  studies  of  unsteady 
flows  at  present. 

6 7 

The  oscillatory  forces  and  moments  found  by  Tijdeman, ' see  Table  1, 
Section  III,  generally  show  less  lag  behind  the  flap  motion  than  the  calculations. 
The  calculated  magnitudes  for  the  forces  and  moments  are  well  scattered 
about  the  experimental  values. 

21 

Kacprzynski,  J.  J. , "Viscous  Effects  in  Transonic  Flow  Past  Airfoils" 

ICAS  Paper  No.  74-19,  The  Ninth  Congress  of  the  International  Council 
of  the  Aeronautical  Sciences,  Haifa,  Israel,  August  25-30,  1974 


Not  much  effort  has  been  expended  on  calculating  flows  where  the  shocks 
have  been  weakened  deliberately  and,  in  controlled  fashion,  shifted  to  positions 
on  the  airfoil  which  agree  with  those  seen  experimentally.  In  steady  transonic 
airfoil  calculations,  palatable-pressure  distributions  might  be  obtainable  by  use 
of  non-conservative  differencing  at  "shocks"  or  by  usii^  coarse  mesh  so  that 
artificial  viscosity  in  the  scheme  distorts  and  smears  the  shock  pressure  rise 
out  over  an  appreciable  fraction  of  the  airfoil  chord.  These  artifices  might  not 
be  useable  for  unsteady  flow  because  they  might  not  properly  model  the  complex 
interplay  between  the  boundary  layer,  the  shock  pressure  rise,  the  shock 
obliquity  near  the  surface,  and  the  shock  propagation  along  the  airfoil  surface 
under  dynamic  conditions. 

Imprecision  in  modelling  the  shock-boundary-layer  effects  might  have  a 
large  effect  on  the  airfoil  pitching  moments  because  the  shock  is  in  the  wrong 
location,  the  motion  of  the  shock  (amplitude  and  wave  form)  during  the  oscilla- 
tion is  incorrect,  and  the  pressure  jump  when  the  shock  crosses  a given  airfoil 
station  is  incorrect.  Unrealistic  shock  motion  is  a feature  of  the  results 
mentioned  in  Section  III  1.4  of  the  present  report.  Definitive  measurements  of 
the  pressure  jump  which  occurs  when  a shock  crosses  a given  pressure  tap  in 
an  experiment  on  an  oscillating  airfoil  would  be  valuable. 

11  14 

Exploratory  calculations  by  Traci  and  by  Ballhaus  and  those  made  in 

the  present  investigation  indicate  a possibility  of  strong  influence  of  the  wind 

timnel  walls  on  results  of  the  quasi  steady  and  unsteady  experiments  of 
6,7 

Ti]  deman. 

The  method  used  here  and  the  others  mentioned  no  doubt  would  be  inade- 
quate if  the  problem  were  to  analyze  transonic  unsteady  flow  over  an  airfoil  with 
partial  or  massive  separation  and  reattachment  occurring  during  the  oscillation. 
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